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Abstract. We use mid-infrared spectral decomposition to separate the 6fim mid-infrared AGN continuum from 
the host emission in the ISO low resolution spectra of 71 active galaxies and compare the results to observed and 
intrinsic 2-lOkeV hard X-ray fluxes from the literature. We find a correlation between mid- infrared luminosity 
and absorption corrected hard X-ray luminosity, but the scatter is about an order of magnitude, significantly 
larger than previously found with smaller statistics. Main contributors to this scatter are likely variations in the 
geometry of absorbing dust, and AGN variability in combination with non-simultaneous observations. There is 
no significant difference between Type 1 and Type 2 objects in the average ratio of mid-infrared and hard X-ray 
emission, a result which is not consistent with the most simple version of a unified scheme in which an optically 
and geometrically thick torus dominates the mid-infrared AGN continuum. Most probably, significant non-torus 
contributions to the AGN mid-IR continuum are masking the expected difference between the two types of AGN. 
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^ 1. Introduction dust, but much more complex radiative transfer effects 

may occur in high optical depth configurations (e.g. Pier 

The mid-infrared and hard X-rays are two of the regions ^ ^ ^^^^^^^ ^^^^^^^ measurements of the mid- 

J3 . of the electromagnetic spectrum that are of particular in- .^^^^^^^ continuum in conjunction with hard X-ray obser- 

■ ■ ' terest for the study of active galactic nuclei (AGN). Hard ^^^-^^^ ^^^^^^^ ^ ^^^^.^^ unification scenarios 

X-rays, unless extremely obscured m fully Compton-thick ^^^^ ^ ^.^^^ ^^j^^.^^ between the two quantities has 

objects, can provide a direct view to the central engine, ^^^^^^j^ ^^^^ reported by Krabbe et al. OOOIJ on the basis 

and are often believed to be a reasonable isotropic mea- ^id-infrared imaging of eight nearby Seyferts. Clavel et 

sure of the bolometric luminosity of the AGN. The nu- ^^^^ ^^^^^ ^ 1^^^^ difference between the equiv- 

clear infrared continuum in Seyferts, in contrast, is due to i .j.i r , i r j itiatt) • • 

' ' alent widths ot the mid-mtrared aromatic FAH emission 

AGN emission reprocessed by dust, either in the putative ^^^^^^^^ ^^p^ ^ 2 Seyferts. Under the assumption 

torus or on somewhat larger scales, e.g. inside the Narrow ^^^^ ^^^^ ^^p^ ^ 2 subsamples are well matched in 

Line Region. The observed mid-infrared emission is thus luminosity and host properties, they interpret this 

a function not only of the AGN luminosity but also of the orientation dependent depression of continuum in 

distribution of the obscuring matter and of the viewing g^y^^^^ 2s with respect to the host-related isotropic PAH 

direction of the observer. In the most simple form this is • ■ 

^ emission. 

due to the covering factor and distance of the obscuring rr^i • • i i i u j.u ^ i ■ i j-£c 

ihe issue is plagued, however, by the technical difh- 

Send offprint requests to: D. Lutz culty of isolating the AGN mid-infrared continuum from 

* Based on observations with ISO, an ESA project with in- t^^e host galaxy emission. IRAS has readily detected large 

struments funded by ESA member states (especially the PI numbers of AGN, but the host contribution to these large 

countries: France, Germany, the Netherlands, and the United beam infrared spectral energy distributions (e.g. Spinoglio 

Kingdom) with the participation of ISAS and NASA. et al. . 1995)1 is not easy to quantify and significant in all but 
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the very powerful AGN. One way to address this difficulty 
is high spatial resolution imaging from groundbased tele- 
scopes, e.g. in the L- and M-bands (e.g. Alonso-Herrero et 
al. lM1T|l or the N-band (e.g. Maiohno et al. lTn^ Krabbe 
et al. l2UUl|l . This method produces reliable results in cases 
of good surface brightness contrast between AGN and 
host, but can face ambiguities in cases where the AGN 
is surrounded by intense star formation on scales similar 
to the spatial resolution used, in particular if the observa- 
tions are diffraction limited by a moderate size telescope 
(e.g. NGC6240, NGC4945; Krabbe et aL EnHTll . 

We use the alternative approach of isolating the AGN 
mid-infrared continuum spectrally, making use of the size- 
able database of low resolution mid-infrared spectra of 
AGN that are a legacy of the Infrared Space Observatory 
ISO. Low resolution spectra of galaxies can be decomposed 
into three components (Laurent et al. l2(J(J(J|l : A component 
dominated by the aromatic 'PAH' features arising in pho- 
todissociation regions or the diffuse interstellar medium of 
the host, an HII region very small grain continuum rising 
steeply towards wavelengths beyond 10/im, and for active 
galaxies a typically flatter thermal AGN dust continuum. 
Starlight is unimportant except for quiescent objects like 
ellipticals or nearby spirals with weak central star for- 
mation or AGN activity. The three components may also 
be obscured, with the additional complication of ice fea- 
tures (Spoon et al. l2(JU2fl . A full spectral decomposition 
accounting for all these effects can be attempted in cases 
of good signal-to-noise ratio and full wavelength cover- 
age (e.g. Tran et al. 120011 Spoon et al. 2003). Since most 
of our data are for the restricted ISOPHOT range (5.8 to 
11.8 ^m) that limits the accuracy of separating silicate ab- 
sorption and PAH emission, and since some spectra are of 
limited S/N, we follow a more straightforward approach. 

In the range covered by the ISOPHOT spectra, the 
AGN emission is most easily isolated shortwards of the 
complex of aromatic emission features (Laurent et al. 
I2000|l . We determine a continuum at 6^m rest wave- 
length, and eliminate non-AGN emission that will in many 
cases be present in the fairly large beam. This is done 
by subtracting a star formation template scaled with the 
strength of the aromatic 'PAH' features arising in the host 
or in circumnuclear star formation. This method does not 
require to spatially resolve the AGN from the contaminat- 
ing star formation, and will face its limits only when trying 
to identify a weak AGN in the presence of strong star for- 
mation or a stellar continuum that can be detectable for 
the most nearby galaxies. Our sample of 71 AGN is then 
used to quantify the relation of mid-infrared and X-ray 
emission at significantly better statistics than previously 
possible. 

2. Data 

We have included in our sample those Seyferts or Quasars 
with low resolution spectra in the ISO archive for which 
hard X-ray observations were available in the literature. 
Since our main goal is a comparison of thermal infrared 



emission with hard X-rays we have excluded objects be- 
lieved to have a significant synchrotron contribution in the 
mid-infrared (e.g. Gen A, 3C273). Of the Ultraluminous 
Infrared Galaxies for which a significant number of ISO 
spectra is available, we have restricted ourselves to a num- 
ber of the brightest objects with clear and undisputed 
AGN contribution. 

2.1. Infrared Spectroscopy 

Most of the spectra are from a rereduction of chopped 
ISOPHOT-S spectra of galactic nuclei in the ISO archive 
(Spoon et al. 2004, in preparation). The reduction was 
done using PIA^ version 9.0.1. Steps in the data reduc- 
tion included: 1) deglitching on ramp level. 2) subdivision 
of ramps in four sections of 32 non destructive read-outs. 
3) ramp fitting to derive signals. 4) masking of bad signals 
by eye-inspection. 5) kappa sigma and min/max clipping 
on remaining signal distribution. 6) determination of av- 
erage signal per chopper plateau. 7) masking or correction 
of bad plateaux by eye-inspection. 8) background subtrac- 
tion by subtracting off-source from on-source plateaux. 9) 
finally, flux calibration, using the signal dependent spec- 
tral response function. The absolute calibration is accurate 
to within 20%. The ISOPHOT-S aperture is 24" x 24" 
and in many cases includes noticeable host emission. A 
few spectra are based on different reductions. Circinus 
and NGC 1068 were observed in staring mode, and for 
NGC 4945 the ramp division of step 2) used only two sec- 
tions. Mrkl014 was reduced with the earlier PIA 7.2 and 
scaled to the PIA 9.0.1 flux scale. For NGC 1808 we have 
used an ISOCAM-CVF spectrum (Laurent et al. l2000l and 
priv. communication). 

We derive the 6/im AGN continuum by a decompo- 
sition over a range from 5.5/xm rest wavelength (or the 
minimum rest wavelength covered) to 6.85/im. With the 
exception of the 6.2/im PAH feature (also associated with 
a continuum- like 'plateau'), this region is relatively free 
of spectral structures or strong emission lines. Absorption 
features from ice (6/im) and hydrocarbons (6.8/xm) are 
found in highly obscured objects but are mostly unde- 
tected in our AGN sample (see Spoon et al. l2002|l . We can 
thus fit the spectrum by the superposition of a star forma- 
tion component dominated by the 6.2/im PAH feature and 
a simple linear approximation for the AGN continuum. 
For the star formation component, we use the SWS spec- 
trum of M82 (Sturm et al. 20D-0j)i rebinned to the exact 
sampling of the individual low resolution spectra. Fig. ^ 
illustrates this using the example of NGC 7469 which is 
well known to host both a powerful Type 1 AGN and a 
strong circumnuclear star forming ring (e.g., Mazzarella 
et al. ll994[ Genzel et al. 11995(1 . The correction for 'contin- 
uum' associated with the aromatic emission is relatively 
modest in this object, but more than half in some Seyferts 
in our sample. The results of the fit procedure are a 6/im 

^ PIA is a joint development by the ESA Astrophysics divi- 
sion and the ISOPHOT consortium 



D. Lutz et al.: Mid-infrared and hard X-ray emission in AGN 



3 




7 8 9 10 11 

Rest wavelength [/j-m] 



12 



O.E 



^ 0.6 



0.4 



0.2 



0.0 



5.8 



6.0 6.2 6.4 

Rest wavelength [^tm] 



6.6 



Fig. 1. Illustration of the decomposition used to isolate 
the AGN continuum. Top panel: ISOPHOT-S spectrum 
of the mixed AGN/starburst galaxy NGC7469, showing 
strong PAH emission as well as an elevated continuum. 
Bottom panel: Cutout of the region around the 6.2/im 
PAH feature. Top continuous line = observed spectrum. 
Top dotted line = fit by the sum of the rebinned SWS 
M82 PAH spectrum and a linear AGN continuum. Bottom 
dotted line = fitted AGN continuum. Bottom continuous 
line — difference of observed spectrum and fitted PAH 
component. The thick vertical line indicates the AGN 6^m 
continuum, obtained by evaluating the fitted continuum 
at that wavelength. An additional correction is applied to 
this value to account for differences between M 82 and an 
average starburst. 



continuum as illustrated in Fig.nand an average flux den- 
sity of the PAH component over the rest wavelength range 
6.1 to 6.35/.tm, both including a 1+z correction so that the 
rest frame spectrum vfi, matches the observed frame one. 

The SWS spectrum of M 82 is technically well suited 
for this decomposition but represents just a single ob- 
ject that may be not fully representative for star form- 
ing objects in general, given a certain dispersion in 
relative importance of 'PDR' and 'HIT components in 
the sense of the model of Laurent et al. H2Ul)U|l . We 
have thus applied the decomposition also to a number 
of bona fide star forming galaxies (NGC23, NGC232, 
NGC 520, IC 342, NGC 3256, M 83, NGC 5653, NGC 6090, 
NGC6946, NGC 7252, NGC 7552), excluding starbursts 



with clear evidence for an additional minor AGN (e.g. 
NGC 3690, Delia Ceca et al. TMV^ . Indeed, M82 is found 
to have relatively low 6/im continuum compared to most 
of these objects. Averaging over all these star forming ob- 
jects including M82 we find a mean residual 6/im contin- 
uum which is 0.096 times the average 6.1 to 6.35/xm PAH 
flux density, with a dispersion of 0.085. We consider both 
this average value and the dispersion in our analysis of 
AGN. The 6/im AGN continua listed in Table [H are the 
values obtained after subtracting from the direct fit result 
0.096 times the average 6.1 to 6.35/im PAH flux density. 
Error estimates for the 6/im continuum are the quadratic 
sum of two components. The first is a measurement error 
based on individual pixel noise derived from the disper- 
sion in the difference of observation and fit. The second 
is 0.085 times the average 6.1 to 6.35/im PAH flux den- 
sity, thus considering the dispersion in the properties of 
the comparison star forming galaxies. Either component 
is found to dominate for part of the AGN sample, de- 
pending on S/N and the relative importance of the PAHs. 
Limits quoted in Tableware 3cr limits based on these con- 
siderations. 

We ascribe the 6/im continuum derived in this way to 
an AGN. This is not strictly correct in all cases, as very in- 
tense star forming environments outside of the parameter 
range covered by our comparison objects also can produce 
a stronger continuum with relatively weak PAH at these 
wavelengths. This is true, for example, for the obscured re- 
gion in SBS 0335-052 (Thuan et al. dQW^ . the mid-infrared 
peak in the Antennae (Mirabel et al. ll998|l . and likely the 
circumnuclear region of Arp 220 (Spoon et al. I2()03|l . We 
consider this of minor importance, however, for our sam- 
ple which does not include star forming dwarfs, and only 
some of the best established AGN among the ULIRGs. 

In some of the weaker AGN in nearby galaxies, pho- 
tospheric emission from the central old stellar population 
may contribute significantly to the measured 6/im con- 
tinuum. The 6/im stellar continumn can be extrapolated 
from the stellar K-band continuum. From four ellipticals 
without PAH emission (NGC 3379, NGC 4374, NGC 2300, 
NGC 4649; Lu et al. 1^71(7^1 Xilouris et aLEOOSI S. Madden 
priv. comm.) we estimate a scaling S'g^ 0.19 X S2. 
We cannot directly apply this extrapolation to our sam- 
ple objects, since it applies only to the stellar K-band, 
and decompositions of the K-band continuum inside our 
aperture into stellar and AGN are usually not available. 
We have verified, however, that the conclusions reached in 
Sect. 3 are robust to stellar continuum contributions. For 
this purpose, we have repeated as an extreme assumption 
our analysis after subtracting directly from the 6/im con- 
tinuum obtained in the spectral decomposition 0.19 times 
the total K-band contimmm in an ISOPHOT-S aperture, 
which can be extrapolated with modest uncertainty from 
2MASS data accessible in NED. The results given in Table 
2 did not change significantly. We did not adopt these val- 
ues, however, since they imply a systematic overcorrection 
for the many AGN dominated objects. Instead, we mark 
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Fig. 2. Distance distribution for the objects in our sam- 
ple, separated by Type 1 (1 to 1.5) and 2 (1.8 to 2) 
Seyferts. Light shading identifies those Seyfert 2's where 
no absorption correction to the AGN hard X-ray emis- 
sion was possible. We adopted distances derived assuming 
Ho=75, qo=0.5 except for M51 (8.4Mpc), M81 (3.63Mpc) 
and Circinus (4Mpc). 



in TableQlthe few cases where a strong stellar contribution 
to our measured 6/im continuum is likely. 

We do not attempt a correction of the mid-infrared 
fluxes for foreground extinction. Published extinction val- 
ues for those objects refer to different tracers, and their 
applicability to the region dominating the mid-infrared 
flux is uncertain. Since the extinction at 6/im is 1/20 or 
less of the visual extinction, typical Ay values from op- 
tical studies are irrelevant in any case. In the discussion 
of individual objects in the Appendix, we also mention 
NGC 6240 and NGC 4945 where obscuration of the mid- 
infrared AGN continuum is likely significant. 

2.2. X-ray emission 

We have compiled the X-ray fiuxes in the 2-lOkeV range 
from various literature sources using data from different 
satellites (e.g. ASCA, BeppoSAX, Chandra, XMM). For 
Seyfert 2s for which the absorbing column density Nh 
could be measured (most cases), the authors generally pro- 
vide reliable measurements of the instrinsic, absorption- 
corrected hard X-ray flux. The intrinsic X-ray flux cannot 
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Fig. 3. Observed hard X-ray fluxes vs. 6/im AGN contin- 
uum for the sample galaxies. Diamonds indicate Type 1 
Seyferts (including objects classified as types 1, 1.2, 1.5), 
asterisks indicate Type 2 Seyferts (including objects clas- 
sified as types 1.8, 1.9, 2). Thin horizontal lines, in most 
cases smaller than the symbols, indicate the la uncer- 
tainty of the 6/im detections. 



be recovered when only a lower limit to Nh is inferred; 
this is the case for totally Compton thick nuclei with 
Nh > lO^^cm"^, which are absorbed at all energies, or for 
mildly Compton thick nuclei with Nh > 10^'*cm~^ which 
are lacking data at E>10keV, i.e. for which any transmit- 
ted flux at high energy cannot be probed. For Seyfert 1 
galaxies without published X-ray column we have assumed 
that observed and intrinsic flux can be safely adopted to be 
equal because of typically low absorbing columns. Figure|21 
shows the distance distributions both for the full sample 
and for the subsample with absorption corrected X-ray 
fiuxes. 



3. Results 

The observed mid-infrared and hard X-ray fiuxes are 
listed in Tableland summarized in Fig.|3 The diagram 
shows the expected separation between Type 1 and Type 
2 Seyferts, hard X-ray emission in the latter often be- 
ing strongly absorbed. More interesting is the comparison 
with absorption corrected X-ray fiuxes shown in Fig. ^ A 
correlation between hard X-rays and infrared continuum 
is indicated in this diagram but with considerable scat- 
ter that is not due to observational error in the mid-IR 
spectra. There is no offset between Type 1 and 2 objects 
any more. Fig. [S] compares mid-infrared and X-ray lumi- 
nosities. Hard X-ray and mid-infrared luminosity correlate 
over four orders of magnitude in luminosity, the ratio of 
the two shows no clear trend with luminosity. The data 
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Table 1. Infrared and hard X-ray fluxes for our sample. The 6/Km flux refers to the AGN continuum exclusively and is 
derived from ISO low resolution spectra as described in the text. Hard X-ray data have been taken from the references 
listed. 



Name 




Type 




S PAH 


JT* 9 1 nt/sl/ nhc 


Nh 


1^9 inti^l/inf 
Z, — L\J tij ff- V 1 1 L t 


Reference 




km/s 




Jy 


Jy 


10~^^ergs^^cm~'^ 


10^"cm"^ 


10~^'^ergs~^cm~^ 




Mrk 334 


6582 


Syl.S 


0.0357 


0.0054 


0.80 


? 




Polletta96 


Mrk335 


7730 


Syl 


0.1508 


<0.0302 


0.96 


4.0 


0.96 


Gondoin02 


NGC 34 


5931 


Sy2 


<0.0804 


0.1810 


<0.39 


? 


? 


Polletta96 


IZwl 


18342 


Syl 


0.1859 


<0.0325 


0.68 


3.0 


0.68 


ReevesOO 


Fairall 9 


14095 


Syl 


0.1386 


<0.0188 


1.53 


2.3 


1.53 


GondoinOla 


NGC 526A 


5725 


Syl.9 


0.1114 


<0.0182 


1.80 


200.0 


1.93 


LandiOl 


Mrk 359 


5215 


Svl 


0.0432 


0.0318 


0.74 


4.8 


0.74 


OBrienOl 


Mrk 1014 


48900 


Svl 


0.0218 


<0.0222 


1.24 


? 


1.24" 


Boller02 


Mrk 590 


7910 


Svl 


0.1132 


0.0232 


1.97 


? 


1.97" 


Malizia99 


NGC 1068 


1148 


Sv2 


11.2195 


1.2491 


0.35 


>le5 




Matt97 


4U0241-I-61 


13200 


Syl 


0.2673 


<0.0269 


3.57 


150.0 


3.82 


Malizia97 


NGC 1097 


1275 


Svl 


<0.3087 


0.8532 


0.17 


6.0 


0.17 


Iyomoto96 


NGC 1144 


8648 


Sy2 


<0.0135 


<0.0119 


<11.00 


? 


? 


PolIetta96 


NGC 1365 


1636 


Syl.8 


<0.5582 


1.2006 


0.66 


4000.0 


2.31 


RisalitiOO 


NGC 1386 


868 


Sy2 


0.1684 


<0.0349 


0.02 


? 




Maiolino98 


Mrk 618 


10658 


Svl 


^0.0179 


<0.0287 


0.70 


3.0 


0.70 


Rao92 


04385-0828 


4527 


Sy2 


0.2205 


<0.0119 


1.80 


? 


?6 


Pollotta96 


NGC 1667 


4547 


Sv2^ 


<0.0542 


0.1360 


0.003 


>le4? 


?5 


Bassani99 


NGC 1808 


1000 


Sv2^ 


<0.5108 


1.9983 


0.09 


320.0 


0.11 


Bassani99 


Ark 120 


9682 


Svl 


0.1106 


0.0215 


3.58 


? 


3.58" 


Nandra94 


05189-2524 


12760 


Sv2 


0.2114 


0.0416 


0.36 


800.0 


0.95 


SevergniniO 1 


MCG8-11-11 


6141 


Syl. 5 


0.1890 


<0.0207 


5.66 


<2.0 


5.66 


PerolaOO 


Mrk 3 


4050 


Sv2 


0.0931 


<0.0177 


0.65 


l.le4 


4.08 


Cappi99 


Mrk 6 


5536 


Syl.5 


0.1424 


<0.0183 


1.20 


1000.0 


2.00 


Immler03 


PG0804-I-761 


30000 


Svl 


0.0778 


<0.0103 


0.88 


<40.0 


0.88 


GeorgeOO 


M81 


-34 


Svl.8 


^0.3835 


<0.0278 


1.50 


9.4 


1.50 


Bassani99 


NGC 3227 


1157 


Svl.5 


0.1518 


0.0862 


0.75 


650.0 


0.76 


Gondoin03a 


UGC6100 


8778 


Sv2 


<0.0186 


0.0249 


<1.14 


? 


? 


Polletta96 


NGC 3516 


2649 


Syl.5 


0.2111 


<0.0635 


4.41 


4.0 


4.41 


GuainazziOl 


NGC 3783 


2550 


Svl 


0.2346 


<0.0210 


8.50 


300.0 


8.50 


Blustin02 


NGC 3982 


1109 


Sv2 


<0.0592 


0.1316 


<0.42 


? 


? 


Polletta96 


NGC 4051 


725 


Svl.5 


0.2014 


<0.0768 


1.60 


<1.0 


1.60 


Reynolds97 


NGC 4151 


995 


Syl.5 


0.9175 


<0.0455 


9.00 


300.0 


9.00 


YangOl 


Mrk 766 


3876 


Syl.5 


0.1415 


0.0354 


3.00 


80.0 


3.00 


MattOO 


NGC 4388 


2524 


Sy2 


0.1298 


0.1213 


1.20 


4300.0 


4.30 


Bassani99 


NGC 4507 


3538 


Sv2 


0.1954 


0.0625 


2.10 


2900.0 


7.03 


Bassani99 


NGC 4579 


1519 


Syl.9 


•^0.0966 


<0.0356 


0.44 


240.0 


0.44 


Bassani99 


NGC 4593 


2698 


Svl 


0.1836 


<0.0271 


3.71 


2.0 


3.79 


Perola02 


NGC 4945 


560 


Sy22 


<0.9092 


2.3230 


0.54 


4.0e4 


19.00 


GuainazziOO 


NGC 5033 


875 


Syl.9 


^0.0691 


0.1588 


0.55 


290.0 


0.55 


Bassani99 


M51 


463 


Sy2 


•^0.0880 


0.1267 


0.03 


5.6e4 


1.28 


FukazawaOl 


NGC 5273 


1089 


Svl. 9 


0.0216 


<0.0190 


<1.17 


? 


? 


Pollctta96 


Mrk 273 


11132 


Sy2 


0.0725 


0.0617 


0.06 


4900.0 


0.35 


Bassani99 


IC4329A 


4813 


Syl 


0.4723 


<0.0250 


16.40 


42.0 


16.80 


GondoinOlb 


Mrk 463 


14904 


Sy2 


0.2518 


<0.0304 


0.04 


1600.0 


0.09 


Bassani99 


Circinus 


439 


Sy2 


4.6688 


4.7199 


1.40 


3.6e4 


40.00 


Matt99 


NGC 5506 


1853 


Syl.9 


0.6607 


0.0665 


8.38 


340.0 


10.80 


Bassam99 


IC4397 


4419 


Sy2 


<0.0228 


0.0587 


2.23 


? 




Polletta96 


NGC 5548 


5149 


Syl.5 


0.1302 


0.0267 


4.30 


51.0 


4.30 


Reynolds97 


NGC 5643 


1199 


Sy2 


0.0610 


0.0926 


0.13 


>lc5 




Maiolino98 


NGC 5674 


7474 


Syl.9 


<n.02r,r, 


0.0562 


0.84 


200.0 


1.28 


Bassani99 


MrkcSli 


i08:)2 


Syi.r) 


O.Ool).") 


O.OilK) 


i.OO 




i.OO' 


l?oyiiokls97 



^ Uncertain classification, sec appendix. 

^ 'Elusive' AGN without NLR detection but listed as type 2 according to the X-ray properties. 
^ Stellar continuum may contribute a large part of the measured 6/um continuum. 
" Assuming negligible absorption for Type 1 objects without Nh measurement. 

^ Intrinsic flux cannot be recovered for Type 2 objects without Nh measurement or with lower limits to Nh- 
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Name czh Type Sg^m Span F2^iokeV,obs Nh F2-iokeV,int Reference 

km/s Jy Jy 10""ergs"^cm"^ 10^°cm"^ 10""ergs"^cm"^ 



NGC 5929 


2492 


Sy2 


^0.0171 


<0.0273 


<0.79 


? 


? 


Polletta96 


15307+3252 


277700 


Sy2 


0.0073 


<0.0023 


<0.007 


? 


? 


Ogasaka97 


PG1613-h658 


38700 


Syl 


0.0394 


<0.0206 


0.69 


3.0 


0.69 


Lawson97 


NGC 6240 


7339 


Sy2 


0.0698 


0.2093 


0.19 


2.2e4 


14.0 


Vignati99 


Mrk 507 


16758 


Sy2 


0.0094 


0.0141 


0.05 


34.0 


0.05 


Bassani99 


H1821+643 


89038 


Syl 


0.0972 


<0.0114 


2.34 


<1.8 


2.34 


ReevesOO 


ES0141-G55 


10793 


Syl 


0.0862 


<0.0238 


2.65 


5.5 


2.65 


Gondoin03b 


19254-7245 


18500 


Sy2 


0.0786 


<0.0593 


0.02 


>le4 


?^ 


Braito03 


Mrk 509 


10312 


Syl 


0.1226 


0.0345 


5.66 


<1.0 


5.66 


PerolaOO 


PKS2048-57 


3402 


Sy2 


0.2581 


0.0153 


1.20 


2400.0 


3.00 


Bassani99 


PG2130-H099 


18630 


Syl 


0.0860 


<0.0227 


0.53 


<19.0 


0.53 


Lawson97 


NGC 7213 


1792 


Syl.5 


^0.1470 


0.0181 


3.66 


? 


3.66'' 


Malizia99 


NGC 7314 


1422 


Syl.9 


0.0442 


<0.0384 


3.56 


120.0 


4.05 


Bassani99 


Ark 564 


7400 


Syl 


0.0652 


0.0134 


2.00 


10.0 


2.00 


TurnerOl 


NGC 7469 


4892 


Syl 


0.1788 


0.3873 


2.90 


<0.4 


2.90 


Reynolds97 


23060+0505 


52200 


Sy2 


0.0963 


<0.0301 


0.15 


900.0 


0.25 


Brandt97 


NGC 7582 


1575 


Sy2 


0.1827 


0.6378 


1.55 


1200.0 


2.72 


Bassani99 


NGC 7674 


8713 


Sy2 


0.1541 


0.0594 


0.05 


>le5 


?^ 


Bassani99 


NGC 7679 


5138 


Syl 


<0.0827 


0.3209 


0.60 


5.0 


0.60 


DellaCecaOl 



o o 



H> NGC1097 

^ NGC 1808 

Mrk463 



0.01 



0.10 1.00 
6|Um AGN continuum [jy] 



10.00 



,46 



,45 



,44 



,43 



,42 



,41 



,40 



,39 



* o * 

NGC624C)-^ 



. O 



M81 NGC 1808 



H1821+643 



O^O 10^1 10^2 ^^43 



AGN continuum luminosity uLi'{6iJ,m) [erg s 



,044 ,o45 ^^46 ^^47 



Fig. 4. Hard X-ray fluxes corrected for absorption vs. 

AGN continuum for those objects where an absorp- 
tion corrected X-ray flux was available. Symbols are as in 
Fig. El The dotted line indicates slope 1, it is not a fit. 



are consistent with the same mean ratio, with considerable 
dispersion, over the entire range of luminosity. 

3.1. The scatter in the relation between mid-infrared 
continuum and hard X-ray emission 

Fig. shows histograms for the ratio of hard X-ray flux 
and mid-infrared AGN continuum, separated by Seyfert 
type. The dispersion is considerable, the ratio for Seyfert 



Fig. 5. Hard X-ray luminosities corrected for absorption 
vs. 6/xm AGN continuum luminosities for those objects 
where an absorption corrected X-ray flux was available. 
Symbols are as in Fig. 13 The dotted line indicates slope 
1, it is not a fit. 

Is varying over an order of magnitude. Ratios for Type 
2 Seyferts scatter yet wider. This may be partly due to 
inaccuracies in the absorption correction in the X-ray 
data. The current data do not allow to determine with 
certainty whether there is in reality a larger spread for 
Type 2 than for Type 1. The two outlying objects Mrk 463 
and NGC 6240 are briefly discussed in the Appendix. The 
spread observed in our sample is significantly larger than 
in the imaging study of Krabbe et al. ((2001,1 . where with 
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E 



Seyfert2 



^3 



-2 -1 

log(F(2-10keV);„,AF„(6Mm)) 




-1 
■10keV);„,AF„(6Mm)) 

Fig. 6. Histograms showing for the two Seyfert types the 
logarithm of the ratio between intrinsic hard X-ray flux 
and 6/xm AGN continuum. Dark shadings indicates ob- 
jects with upper hmits for the 6/im continuum, i.e. lower 
limits for the ratio shown. 



the exception of the outlier NGC 6240 the ratio of intrinsic 
hard-X and mid-IR emission varies by just a factor of 3-4. 
This was probably a fortuitious effect of the small sample 
including ~10 times less objects than the present study. 
Given the many factors of AGN spectral energy distribu- 
tion and geometry relating the hard X-ray flux and the 
dust reradiation of AGN emission in the mid-infrared, the 
larger scatter is not surprising. Another signiflcant contri- 
bution to the spread is likely due to AGN variability. The 
X-ray and infrared measurements summarized in Tabled 
were usually taken several years apart. Even in case of 
simultaneous observations, short term variability of the 
central engine may contribute to the scatter because of 
different time response and time averaging effects for the 
X-rays and for the larger scale dust. 

3.2. No significant difference between Seyfert 1 and 
Seyfert 2 types 

A surprising result given the reasonable statistics of the 
present sample is the failure to detect the difference be- 
tween Type 1 and Type 2 Seyferts that is predicted by 
the most straightforward versions of unified models. The 
median log(F(2-10keV)/i/F^(6/im)) is -0.41 for Type 1 ob- 



jects and -0.63 for Type 2. This would not change strongly 
replacing the few lower limits by detections (Fig. ^ . The 
key point is not the small difference found, which might 
even be reduced slighlty if some of the Sy2 limits were 
replaced by detections or with full correction for stellar 
continuum, it is the failure to detect a strong and sig- 
nificant difference in the opposite direction. If the mid- 
infrared continuum in Seyfert 2s is suppressed by a factor 
of ~8 (Clavel et al. 1210111), the hard X/IR ratio should be 
higher by the same factor compared to Seyfert Is. While 
our targets are from a heterogeneous set of ISO observing 
programs, they do not represent a preselection by mid-IR 
flux which may affect such a comparison. The main ob- 
serving programs involved did not invoke such a selection. 
A major part, e.g., is formed by objects from the CfA 
sample (Huchra & Burg 1199^ . Also, for example, most 
AGN from the hard X-ray selected sample of Piccinotti et 
al. ( 19821 are included. As discussed by Maiolino & Rieke 
(,1995^1 . samples like the CfA one may be biased against 
obscured objects and thus not reproduce the real fractions 
of Seyfert types. Our analysis normalizing to intrinsic X- 
rays should be robust to such effects as long as reaching 
lower but still signiflcant numbers of obscured systems. 

In unified schemes (e.g. Antonucci et al. I1993|l . the 
difference between Seyfert types is due to effects of view- 
ing intrinsically similar objects from different directions, 
because an anisotropic distribution of absorbing material 
(e.g. the 'torus') absorbs, scatters, and reprocesses the 
direct AGN light. In the most simple version, a central 
very small source (also emitting the X-rays) illuminates 
a torus-like geometrically and optically thick dusty struc- 
ture. Radiative transfer models of the emission from a 
geometrically and optically thick torus (e.g. Pier et al. 
119921 Efstathiou & Rowan-Robinson 119951 Granato et al. 
119971 Nenkova et al. I2[)()2(l predict a strong anisotropy 
of the mid-infrared emission. The unified model in this 
case clearly predicts for Seyfert 2 galaxies a higher ratio 
of absorption corrected X-ray emission and mid-infrared 
continuum. 

In a simple unified scheme invoking a torus, the ratio 
of intrinsic hard X-ray to observed mid-IR emission will 
increase from a fully face-on to a fully edge-on view. The 
difference between Type Is and Type 2s could hence be 
masked if the Type 2s in the sample were preferentially 
objects with still relatively low obscuring columns, likely 
seen at intermediate angles. While we had to exclude fully 
Compton-thick Sy2s from the comparison because of the 
impossibility to derive an absorption corrected X-ray flux, 
they form well below half of the Sy2 sample (5 of 38, more 
could be among the 11 Type 2 objects where the absorbing 
column is undetermined or only a limit to the 2-lOkeV flux 
available). Abandoning for a moment the normalization to 
X-rays, we have compared the average ratio of 6/im AGN 
continuum to 6.2/im PAH feature for the Seyfert 2's with- 
out intrinsic X-ray fiux in our sample and the ones where 
an intrinsic fiux could be derived. Both using the measure- 
ments/limits in Tableland using the raw measurements, 
the median AGN continuum to PAH ratio is only insignif- 
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NGC6240 



1 



-,20 



-,21 



,22 



,23 



,24 



25 



10^"" 10^ 10^^ 10^'^ 10^ 10 

X — ray absorbing calumn [cm~^] 

Fig. 7. Ratio of intrinsic hard X-ray flux and 6/im AGN 
continuum, plotted as a function of X-ray absorbing col- 
umn. 



icantly (^20%) smaller for the 16 Type 2 objects without 
intrinsic X-ray flux. With all caveats of normalizing to the 
PAHs rather than to intrinsic X-rays, this tentatively in- 
dicates that in a complete Sy2 sample including the fully 
Compton-thick objects the mid-infrared AGN continuum 
will not be much lower. 

More insight can be obtained from plotting the ratio of 
X-ray and mid-infrared emission versus the X-ray absorb- 
ing column (Fig. Uj). Column densities above lO^'^cm"^ 
are well populated but do not show the upturn in the 
hard X/Mid-IR ratio expected at high columns in the sim- 
ple unified scheme. In a naive foreground screen the op- 
tical depth at 6/im exceeds one already below lO^^cm^^. 
Radiative transfer in different more realistic 'torus' geome- 
tries (see the large range in the models cited above) will 
predict different magnitudes of this effect and place it at 
different column densities. For example, a strong and fairly 
sharp upturn at Nh 6 x lO^'^cm"^ is expected for the 
best fitting tapered disk model 5e of Efstathiou & Rowan- 
Robinson flOOSI. Placing predictions of future torus mod- 
els on a fully empirical diagram like Fig. [3 should provide 
a valuable test of such models that does not make direct 
assumptions on viewing angles. Observationally, popula- 
tion of the columns around and above 10^^cm~^ will be 
important to characterize the effect of missing the most 
absorbed objects in a sample like ours. 

One possibility why we fail to observe the expected 
signature of a torus is that much of the mid-infrared emis- 
sion is from a more extended region emitting isotropi- 
cally, and overwhelming a possible anisotropic emission. 
Observations have shown this to be the case for the lO/im 
emission of NGC 1068 where at least two thirds of the 
emission comes from an extended region overlapping the 



0) 




38 40 42 44 

bg(6.2/xm PAH luminosity [erg s~']) 



E 4 



Seyfertl 



38 40 42 44 

og(6.2/.im PAH luminosity) [erg s~'] 

Fig. 8. Distribution of luminosities of the 6.2/zm PAH 
emission feature. Dark shading indicates upper limits. 



Narrow Line Region (e.g., Cameron et al. ll99HI Bock et al. 
1998) rather than a compact parsec-scale torus. A similar 
result (at least 27% extended emission) has recently been 
found for NGC 4151 (Radomski et al. IMlHjl . High reso- 
lution observations and mid-infrared interferometry with 
large telescopes should be able to test in the near future 
whether the same applies to other Seyferts. This scenario 
of (mostly) larger scale and isotropic AGN continuum at 
first glance contradicts the results of Clavel et al. (20001), 
however, who find in their large aperture ISOPHOT-S 
data a ^^8 times lower AGN continuum to host PAH ratio 
in Seyfert 2s compared to Seyfert Is. They ascribed this to 
orientation effects in a unified scenario, of a magnitude sig- 
nificantly larger than suggested earlier by Heckman H1995|l 
and Maiolino et al. (.1995)) . Much of their result could be 
an AGN luminosity effect rather than an orientation ef- 
fect, though. Some of the relevant diagnostics are shown 
in Figs. 151 to 1101 and Table [21 for our sample which has a 
significant overlap with the Clavel et al. H2000() sample. 
The distributions of PAH luminosities for the two types 
are consistent (Fig. |Hl, but the caveat that only upper 
limits to the 6.2/im PAH are available for many of the 
strong continuum Syls has to be noted. This result means 
there are no obvious differences in the host properties as 
probed by the PAH emission, although some differences 
(as suggested e.g. by Maiolino ct al. 11995(1 might be hidden 
by the many PAH upper limits, or by sample and phys- 
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Table 2. Properties of the Type 1 and Type 2 objects. The numbers quoted refer to the combination of detections 
and hmits (see Figs. |Hlto^|for the fraction of Hmits which is highest among the Sy 1 PAH luminosities). Distances 
refer to the full sample, all others to the part with extinction corrected hard X-ray data (32 Type 1, 23 Type 2) 



Property 


Seyfert 1 


Seyfert 2 




Median 


Disp. 


Median 


Disp. 


log(Distance) [Mpc] 


2.00 


0.50 


1.68 


0.61 


log(6.2/im PAH luminosity) [erg s~^] 


(41.88) 


(0.85) 


41.50 


0.99 


log(AGN continuum luminosity i^L;^(6/im)) [erg s^^] 


43.88 


0.84 


43.20 


1.00 


log(Extinction corrected 2-lOkeV luminosity) [erg s~^] 


43.39 


0.99 


42.64 


1.07 


log(F(2-10keV) //xF^ (6^m) ) 


-0.41 


0.40 


-0.63 


0.69 



E 4 
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log(Absorption corrected 2-10keV luminosity) [erg s~'] 
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Seyferti 





40 42 44 46 
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Fig. 9. Distribution of luminosities of the 6/im AGN con- Fig. 10. Distribution of extinction corrected hard X-ray 



tinuum. Dark shading indicates upper limits. 



luminosities. 



ical aperture size (distance) mismatch between the two 
types. We note that upper limits also apply to the 6.2 and 
7.7/im PAH features in many objects presented by Clavel 
et al. I|2()()()|l . rather than detections reported through their 
method of including only positive flux above the noisy con- 
tinuum into their PAH fluxes. For our sample, both the 
mid-infrared continuum luminosities (Fig. [SJ and the ab- 
sorption corrected X-ray luminosities (Fig. I1U|I are larger 
in Seyfert Is than in Seyfert 2s. Reflecting the luminosity 
differences, a difference is also seen in the distance dis- 
tributions (Fig. 12)). All that argues that a difference in 
PAH equivalent width between Seyfert Is and 2s, found 
by Clavel et al. (2000), and also seen in our data (com- 
pare Fig. inland Fig.EJ, can be largely due to the Seyfert 
1 part of the samples representing both higher AGN lu- 



minosities and higher AGN/host luminosity ratios. This 
is a natural consequence of many AGN samples, includ- 
ing ones driving the selection of ISO observations, being 
biased in favour of more distant Type 1 objects (see discus- 
sion in Maiolino & Rieke ITM^ . Explaining our results by 
a dominant component of isotropic mid-infrared emission 
thus remains viable. We emphasise that a minor strongly 
anisotropic mid-infrared component may be present with- 
out affecting our analysis. 

Another factor masking the signature of unification in 
the mid-IR to hard X-ray ratios could be uncertainties in 
the absorption correction of the X-ray fluxes. While we 
have argued above that they may contribute to the in- 
creased scatter observed for the Type 2 objects, we see 
no evidence for a systematic offset, which would also have 
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profound implications on estimates of the energy budget of 
obscured X-ray sources. Similar arguments apply to scat- 
ter induced by intrinsic variations of the AGN spectral en- 
ergy distributions between the hard X-rays and the wider 
spectral range that is heating dust. Finally, the strategy 
of our comparison relies on the assumption of isotropic 
X-ray emission (once considering the absorption correc- 
tions). If X-rays were emitted anisotropically and if there 
were alignment between the structures causing this X-ray 
anisotropy and a torus absorbing/emitting anisotropically 
in the mid-IR, the expected difference in the ratio of hard 
X and mid-IR emission could cancel largely. 

Our results suggest that local, mostly moderate lu- 
minosity, Seyferts do not show the behaviour expected if 
their mid- infrared emission were dominated by a compact, 
anisotropically emitting torus. This adds to evidence that 
the most simple version of a unified scheme, in which all 
AGN at all luminosities are governed by such a structure, 
is not applicable. Recent deep X-ray and mid-IR surveys 
indicate significant differences in the evolution of high lu- 
minosity Type 1 and lower luminosity Type 2 populations 
that cannot be reconciled with a single type of obscuring 
structure in AGN of all luminosities and redshifts (e.g. 
Franceschini et al. l2(JU2fl . 

4. Conclusions 

We have used spectral decomposition of a large sample 
of ISO spectra of AGNs to isolate the AGN 6/im con- 
tinua. We compare these mid-infrared continua to intrin- 
sic hard X-ray fluxes from the literature, assumed to be a 
fair isotropic measure of AGN luminosity. Due to this nor- 
malization, our comparison can test AGN properties and 
unification aspects without the level of sensitivity to selec- 
tion biases (e.g. luminosities, ratio of Type 1 and Type 2 
objects in the sample) that is found in comparisons of ab- 
solute quantitities or of normalizations to non-AGN quan- 
titities. The main results are: 

(1) Mid-infrared and intrinsic X-ray fluxes correlate, 
but with signiflcant scatter. Ratios vary by more than an 
order of magnitude, the dispersion in the log is 0.4 for 
Seyfert 1 and 0.69 for Seyfert 2 galaxies. Main contributors 
to this spread likely include variations in the AGN spectral 
energy distribution and geometry of obscuring dust, as 
well as the effects of AGN variability. 

(2) There is no signiflcant difference between Type 
2 and Type 1 in the average ratio of X-ray and mid- 
infrared continuum, in contrast to expectations from uni- 
fied scenarios invoking an optically and geometrically 
thick torus emitting anistropically in the mid-IR. Most 
likely, this is due to a large contribution from extended 
dust components emitting more isotropically. This will 
dilute anisotropic emission that, however, may still be 
present at lower levels. 
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Appendix A: Notes on selected objects 

NGC 1808 is believed to host both a strong spatially ex- 
tended starburst (e.g. Krabbe et al. I1994|l and possibly 
a weak and likely fading AGN (e.g. Bassani et al. I1999|l . 
although the X-ray luminosity is low enough to overlap 
the regime of (also variable) non-AGN Ultraluminous X- 
ray Sources (ULX). Such objects are found at non-nuclear 
positions in other starbursting objects (Fabbiano et al. 
I2003|l . Because of the combination of strong star forma- 
tion and weak AGN, the limit on AGN continuum is far 
off the correlation between mid-IR and X-rays. 
NGC 1667 is likely another example of a fading 'fossil' 
AGN (Bassani et al. I1999|l , introducing additional uncer- 
tainty on the AGN contribution at different wavelengths. 
For example, if the AGN is turning off the apparently high 
Nh may be an artifact of the reflected component remain- 
ing visible longer. 

NGC 4945. The 6/im continuum is high compared to the 
average starburst, but not at the level of a signiflcant de- 
tection of the AGN continuum given the scatter in star- 
burst properties. The mid-IR weakness and correspond- 
ing location in our diagrams is certainly in part due to 
the foreground extinction which is high even in the mid- 
infrared in this system (Spoon et al. JOOOj. A high reso- 
lution mid-IR image, if possible outside the silicate band, 
is needed to break the ambiguity of both this result and 
the low resolution mid-IR imaging of Krabbe et al. l|2001|l . 
and reach a mid-IR detection of the AGN. 
NGC 6240 is one of the sources deviating most from the 
correlation, in the sense of high X-ray and low mid-IR. 
This could include a contribution of adopting a high es- 
timate for the intrinsic X-ray emission (cf. also the lower 
estimate of Ikebe et al. l2000|l . and the likely noticeable ob- 
scuration of the 6/im continuum in this object. See Lutz 
et al. H2003|l for a full discussion of the mid-infrared prop- 
erties. 

Mrk463 has a very strong mid-IR AGN continuum but 
the relatively lowest absorption corrected hard X-ray emis- 
sion compared to the average relation. The X-ray column 
may be underestimated. 
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